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ABSTRACT 

Ultravoilet (UV) absorption lines provide abundant spectroscopic information enabling the probe of 
the physical conditions in AGN outflows, but the outflow radii (and the energetics consequently) can 
only be determined indirectly. We present the first direct test of these determinations using integral 
held unit (IFU) spectroscopy. We have conducted Gemini IFU mapping of the ionized gas nebulae 
surrounding two AGNs, whose outflow radii have been constrained by UV absorption line analyses. 
In Mrk 509, we And a quasi-spherical outflow with a radius of 1.2 kpc and a velocity of ^ 290 km s“^, 
while IRAS F04250—5718 is driving a biconical outflow extending out to 2.9 kpc, with a velocity of 
^ 580 km s“^ and an opening angle of ^ 70°. The derived mass flow rate is ^ 5 and > 1 Mq yr“^, 
respectively, and the kinetic luminosity is > 1 x 10^^ erg s“^ for both. Adopting the outflow radii and 
geometric parameters measured from IFU, absorption line analyses would yield mass flow rates and 
kinetic luminosities in agreement with the above results within a factor of ^ 2. We conclude that the 
spatial locations, kinematics and energetics revealed by this IFU emission-line study are consistent 
with pre-existing UV absorption line analyses, providing a long-awaited direct confirmation of the 
latter as an effective approach for characterizing outflow properties. 

Subject headings: quasars: individual (Mrk 509, IRAS F04250—5718) — quasars: emission lines — 
quasars: absorption lines 


1. INTRODUCTION 

The active phase of the super-massive black holes re¬ 
siding in galaxy centers has crucial effects on the global 
properties of their host galaxies. Thus, feedback from 
active galactic nuclei (AGNs) has become an indispens¬ 
able ingredient of galaxy formation models, requiring the 
black holes to blow large-scale outflowing winds which 
clear out the s urrounding gas and prevent future star 
form a tion fe.g. [ Tabor fc BinnevI 119931: iDi Matteo et al.l 
l2nM IHonkins fc ElvifJ 1201 OH . Roughly 20-40% of 
quasars show blueshifted broad absorption lines (BAL), 
implyin g that snb-relatiyistic q u asar outflows are ubiq¬ 
uitou s iHewett fc Foltzl I2QQ3I: iGangu R fc BrothertonI 
I2QQ8I: iKnigge et al.l I2QQ8I: iDai et al.l l2QQ^ . Theo¬ 
retical modeling has shown that massive, wide an¬ 
gle outflows can indeed be ver y efficient fe e dback 
agents fe.g. I Ostriker et al l I2Q1QI: iCiotti et al.l 2010 : 


t s (e.g. lUstriker et al I I2D1UI: KJiotti et al.l ; 

r]l2QlQUMcCarthv et al.l2Ql(l Hopkins fc Elvii2QlQ: 


Soker 


Faucher-Gignere et al.l 1201^ iChoi et aLll2Q14f). 

Over the last decade, our collaboration and other 
research groups have been performing absorption 
line analyses on about 20 AGN absorption out¬ 
flows to measure the distance from the central 
sour ce, the mass flow rate and kinetic lumino s¬ 
ity (iHamann et ^ I2QQ1I: Ide Kool et al.l l2QQlL l2QQ2allic 


Gabel et al.ll2QQ5l: iMoe et al.ll2QQ9l: iBantista et al.l 

Dunn et al ] i2nRit [yn et al 1 HquIJa rav et all 


Borguet et al.l l2Q12allbl. I2Q13I: lEdmonds et al.l 


Arav et al.l I2Q13I: lLu(w et ^ l2Ql4 iChamberlain et af] 

2Q15[) . Our state-of-the-art analysis methods are able to 
determine the outflow sizes to better than a factor of two 
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(jBautista et al.ll2QlQF frequent ly finding them located at 
kpc scales (see lArav et al.ll2QT^ for a review), comparable 
to the spatial extent of the host galaxies. The main weak¬ 
ness of these investigations, obtained by ultraviolet (UV) 
spectroscopy of excited transitions coupled with detailed 
photoionization modeling, is the lack of spatial informa¬ 
tion. The indirectly inferred outflow sizes remain to be 
confirmed, not to mention their morphology and spatial 
structure. 

In recent years, integral held unit (IFU) observations 


have been c onducted to map the i onized gas surrounding 


/AUriNs le.g. iiNesvaaoa ex ai 

Fu & StocktonI 120091: iLioai 

.Mzuuoi: inusemann ex ai.iizuuoi: 

4 et al.l 120091: 1 Alexander et al. 

2010|: 

iHumohrev et al.l 

[20T( 

It iRuoke & Veilleux 

2^: 

iHusemann et al.l 

I 2 OIII: ICano-Dfaz et al. 


iHarrison et al. 

120121: 

iWestmoQuette et al. 


IHusemann et al.l 

ioil 

iLiu et al.l l2013allbl: 

Ruoke 

& VeilleuxI 120131: It 

nrrison et al.l 120141: iLiu et al.l 

2014 

iShih & StocktonI 

20141; 

iMcElrov et al.l 120151: 


enabled us to detect massive, powerful, kpc-scale, 
AGN-driven outflows both in ra dio-quiet major mergers 
(jRupke fc Veilleu^ l2QllL 1201311 an d in fl igfl luminos¬ 
ity radio-qui e t typ e 2 (|Liu et al.l l2Q13aro and type 
1 (jLiu et al.l l2Q14l ) quasars, providing evidence for 
ubiquitous galaxy-wide feedback in luminous quasars. 

IFU spectroscopy has the obvious advantage of directly 
measuring the spatial extent of the outflows which can 
only be derived indirectly from absorption line analy¬ 
ses (not to mention the morphology and the spatial dis¬ 
tribution of gas kinematics that can only be obtained 
with IFU), yet the latter provides more abundant spec¬ 
troscopic information covering a wide range of ioniza¬ 
tion and density, enabling the probe of the physical con- 
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dition s in the multi-phase outflowing gas (|Arav et al.l 
l2Q13f ). Therefore, IFU mapping of outflowing AGNs will 
directly test the outflow locations derived from UV ab¬ 
sorption lines, providing a strong constraint on the phys¬ 
ical assumptions and parameters of the outflow models. 

In this paper, we present Gemini IFU mapping of the 
ionized gas nebulae around two bright radio-quiet ac¬ 
tive galaxies, Mrk 509 and IRAS F04250—5718. The 
locations of their outflows both have been determined 
or constrained by previous UV absorption line analyses 
using high-quality spectroscopy. 

Although spectroscopically classified as Seyfert 1.5 
galaxies (Table [T]), Mrk 509 and IRAS F04250—5718 
both reside close to the q uasar/Seyfert galaxy border - 
line (|Kopvlov et al.l 11974 IVeron-Gettv fc VeronI l2006f ) , 
and are thus deemed low-luminosity quasars. As one 
of the best studie d AGNs, Mrk 509 {z = 0.03440 
iHuchra et al.l 1 19931) has been scrutinized in multiple 
windows of its electromagnetic spectrum (a thorough 
de scription of its obse rvational properties is available 
in iKaastra et al.l l201l[ ). Combining HST/COS data 
fro m a massive multi-wa velength monitoring campaign 
fsee IKaastra et al1l201l[ ) and archival HST/STIS data, 
lArav et ah (j2012[ ) analyzed all the kinematic components 
of its outflow and obtained a lower limit of 100-200 pc 

for their galactocentric distance. __ 

IRAS F04250-5718 {z = 0.104. PTh^as et al.lIT^ 
was discovered in X-ray as IH 0419—577 IWood et al.l 
11984 and is also known as LB 1727 and lES 0425—573. 
So far, the X-ray spectral energy distribution has 
been the primary interest of the existing studies on it 
(the reader is r eferre d to the introduction section of 
lEdmonds et al ] (I 2 M 1 f or a detailed de s criptio n of its 
observational history). lEdmonds et al.l (j2011f ) identi¬ 
fied three kinematic components from its HST/COS UV 
spectra, and determined a conservative lower limit of 3 
kpc for the outflow. The basic characteristics of these 
two objects and the observing information of our obser¬ 
vations are summarized in Table [H 

This paper is organized as follows. In Section [21 we de¬ 
scribe the observations and reduction of our data. Elux 
calibration, PSE subtraction and line fitting are pre¬ 
sented in Section [3l Outflow features are identified and 
characterized in Section IH In Section [5l we measure 
the spatial scales and derive physical velocities of the 
outflows, and compare to absorption line analyses. The 
mass flow rates and kinetic luminosities are estimated 
in Section [6l followed by a summary in Section [71 We 
adopt an h = 0.70, = 0.29, Oa = 0.71 cosmology 

throughout this paper. 


2. OBSERVATIONS AND DATA REDUCTION 

We observed the program targets with the Gemini- 
South telescope between 2013 July and 2014 January 
(program ID: GS-2013B-Q-84, PI: D. Rupke; we also ob¬ 
served IRAS E22456—5125 in our campaign, but the data 
quality is insufficient for deriving useful results). We used 
the 1-slit Integral Eield Unit (lEU) mode on GMOS-S, 
so that the wavelength coverage is optimized. We per¬ 
formed two exposures, separated by a spatial dithering 
of 1.5" along the shorter axis of the field of view (EoV). 
Combining the dithered exposures, we cover a final EoV 


of 5" x5", translating to a linear scale of 3.5x3.5 kpc^ for 
Mrk 509 and 9.5 x9.5 kpc^ for IRAS E04250-5718. The 
science EoV, sampled by 1000 contiguous 0.2"-diameter 
hexagonal lenslets, is interpolated to a finer grid (0.1") 
in the final data products. 

The seeing at the observing site was 0.6"-0.7", mea¬ 
sured from the point spread function (PSE) model we 
construct from the continuum emission (see Section [3^ . 
We used the spectroscopy blocking filter GG455-G0329 
in our program. The central wavelength of the em¬ 
ployed grating B600-G5323 was set to 625 nm for IRAS 
E04250—5718 and 615 nm for Mrk 509, respectively. This 
grating has a spectral resolution of R ^ 2300 at these 
wavelengths. Both objects were observed between 4700A 
and 7600A so as to cover the [O lll]-H/3 region. To ac¬ 
curately determine the spectral resolution, we perform 
Gaussian fits on unresolved sky lines in the [O lll]-H/3 
neighborhood, finding their full width at half maximum 
(EWHM) to be ^ 120 km s“^. The emission lines of 
interest ([O iii]AA5007, 4959A, H/3, Ha, [N ii]AA6548, 
6583A) are all well resolved. Table [T] summarizes our 
Gemini observations. 

We use the Gemini package for IRAF0 to reduce the 
data, following the standard procedure for CMOS lEU 
described in the tasks gmosinfoifu a nd gmosexamjjle s 
with slight modification (detailed in iLiu et al.l l201^ . 
The datacube product reduced from each exposure has 
0.1" spatial pixels (“spaxels”). We finally combine the 
two frames using the task imcombine by taking the mean 
spectra in each spaxel. The [O lll]5007A line is well de¬ 
tected in every spatial position, with signal-to-noise ra¬ 
tios at its peak of 10-20 at the EoV edges, and 2500-2800 
in the EoV centers. 

3. IFU DATA ANALYSIS 

3.1. Flux calibration and continuum / Fe ll subtraction 

We flux-calibrate our data against previous spec- 
troscopic measuremen ts. Eor IRAS E04250—5718, 
iGuainazzi et al.l (j 19981 ) find an integrated [O ill] flux of 
/[Oiii] = 5.0 X 10“^^ erg s“^ cm“^ using the ESO 1.52m 
telescope equipped with an 8"-wide long slit, while a 5- 
year coordinated program of spectroscopic monitoring of 
Mrk 509 reports /[om] = 0.79 x 10“^^ erg s~^ cm“^ using 
various apertures si zes, of which the smallest is 5" x 8" 
(jCarone et al.lll99^ . Elux loss from these measurements 
are deemed negligible because the employed apertures 
are larger than our field of view. 

In order to remove the contamination from Ee II emis¬ 
sion, we perform a procedure identical to iLiu et aP 
(|2014j ). We fit the continua in the [O lll]-H/3 neighbor¬ 
hood with the sum of two components: a quadratic poly¬ 
nomia l, and the Ee ll template from iBoroson fc Gree^ 
(jl992f ) smoothed using a Gaussian kernel whose width is 
one of the free fitting parameters. However, we note that 
Ee II contamination to our [O ill] analysis is insignificant 
in both targets. 

3.2. PSF subtraction 

^ The Image Reduction and Analysis Facility (IRAF) is dis¬ 
tributed by the National Optical Astronomy Observatories which is 
operated by the Association of Universities for Research in Astron¬ 
omy, Inc. under cooperative agreement with the National Science 
Foundation. 
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Table 1 

Characteristics of the sample objects and observational parameters. 


Object 

a, S (J2000.0) 

z 

Type 

PA 

texp 

Seeing 

L[o III] 

^[O ml, peak 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

Mrk 509 

20‘’44“09?7 -10°43'24"5 

0.0344 

Seyfert 1.5 

0° 

1200x2 

0{'70 

1.9x10^2 

4.9x10-13 

IRAS F04250-57I8 

04''26“0a:7 -57°12'01''0 

0.104 

Seyfert 1.5 

90° 

1600x2 

0''62 

1.4 

2.5x10-13 


Note. — (1) Object name. (2) Right asc ension and declination. (3) Redshifts are taken from lThomas et al ] (ESHl, for 
IRAS F04250—5718) and IHnchra et ahl (|1993l . for Mrk 509). (4) AGN type, from the NASA/IPAC Extragalactic Database 
(NED, http://ned.ipac.caltech.edu). (5) Position angle of the field of view, measured by the longer (5^^ of FoV 
(1-slit 14'U) from north to east. (6) Exposure time (in seconds) and number of exposures. (7) Seeing at the observing site 
(FWHM). (8) Total luminosity of the [O lll]A5007 A line (erg s“^), from observations bv IGuainazzi et ^ (119981 . for IRAS 
F04250—5718) and I Car one et ahl (|1996l . for Mrk 509). (9) Peak [O lll]A5007 A surface brightness before PSF subtraction 
(erg s“^ cm“^ arcsec“^). 


Our program was undertaken when seeing at the ob¬ 
serving site was 0.6"-0.7", requiring careful analysis of 
our datasets to minimize the effects of the PSF. The ob¬ 
servations become significantly easier to interpret if the 
PSF is carefully modeled and subtracted from the IFU 
datacube. The most reliable approach is constructing 
PSF models from the real data themselves. 

Quasar light sca ttered by the interstellar matter 
( Borguet et al.[2QQ^ or star formation in th e quasar host 
( Letawe et al.l 20071 : iSilverman et al.l l2009[) might both 
contribute to the continuum emission. However, we as¬ 
sume that the continuum emission is dominated by a 
central point source in this work. This is a reasonable 
assumption for both of our targets that are classified as 
type 1 AGNs, and is further validated by the fact that 
the FoV is only larger than the FWHM of the PSF by 
a factor of 6, and thus the scattered light from the host 
galaxy, if detectable at all in the outer regions, is over¬ 
whelmed by the broad wings of the PSF. 

To construct a PSF model in an optimized way, we 
first select three rest-frame wavelength intervals free of 


[O III] and H/3 line emission, 4925-4945A, 4970-4980A 
and 5020-5050A, normalize the central brightest pixel 
to unity, then create a median image for each interval. 
For where [O ill] and H/3 emission exist, we linearly in¬ 
terpolate between each pair of adjacent median images 
to construct the PSF model at each wavelength, so that 
a PSF model datacube is created. These constructed 
PSFs have an FWHM of 0.62" for IRAS F04250-5718 
and 0.70" for Mrk 509. 

For the [O lll]-H/3 region, we perform Gaussian fits 
to both the PSF model at each wavelength and to the 
corresponding slice image to compute their centroids us¬ 
ing the GGNTRD routine in the IDL Astronomy User’s 
Librar3l3 (Gaussian fits are sufficient for our purpose of 
centroiding only). We then align the PSF to the slice im¬ 
age using the obtained centroids (although the position 
of the AGN is fixed in the sky, we allow the peak of the 
PSF to move minimally in the plane of the sky (within 
one spatial pixel, 0.1") for optimized results). After that, 
the PSF is scaled to the central peak brightness of the 
slice image, and subtracted from the latter. As the last 
step, guided by these automated algorithms, we manu¬ 
ally adjust the brightness of the PSF to further correct 
imperfect subtractions. 


^ http://idlastro.gsfc.nasa.gov/ 


3.3. Spectral line fits 

After removing Fe ll emission and subtracting the PSF 
from the IFU datacubes, we perform multi-Gaussian fits 
to the [O III] 5007,4959 doublet, so that a noiseless model 
of the line i n every spatial pix el is attained, a procedure 
identical to iLiu et al.l (j2013b[ ). As described in that pa¬ 
per, up to 3 Gaussian components are needed for good 
fits so that the reduced < 2 in every spatial pixel, 
and the actual number of employed Gaussians is deter¬ 
mined by comparing their respective reduced values. 
We then compute the [O ill] line intensity in every spa¬ 
tial position from the multi-Gaussian fit (instead of the 
observed profile). The surface brightness sensitivity of 
our [O III] surface brightness maps is measured to be 
^[o III] = (3-6) X 10“^^ erg s~^ cm“^ arcsec“^. 

As in iLiu et all (j2Q13b[ ). we do not assign physical 
meanings to the individual Gaussians (with the sole ex¬ 
ception of decomposing [O ill] to r emove the tidal tail 
emission in Mrk 509, see Section|4]2j). The fits are merely 
used to produce noiseless models of the [O ill] line to fa¬ 
cilitates non -parametric me asure ments. 

Following IWhittld (jl985[ ) and ILiu et al] (j2Q13b[ ) , we 
measure two quantities that characterize the line-of-sight 
velocity and velocity dispersion of the ionized gas, the 
median velocity (umed) that bisects the total area under¬ 
neath the [O III] emission line profile, and the velocity 
interval that encloses 80% of the total [O ill] emission 
centered at the median velocity (Hso)- As noted in those 
papers, Wgo is more sensitive to the weak broad wings 
of a non-Gaussian profile, but is similar to FWHM for a 
Gaussian velocity profile {Wso = 1.088 x FWHM). 

4. OUTFLOW SIGNATURES 
4.1. IRAS F04250-5718 

The [O III] surface brightness, median velocity and 
Wso maps are shown in Figure [TJ A pair of strong 
[O III] -emitting blobs reside in both sides of the central 
AGN (the central black pixels are due to slight over¬ 
subtraction), featuring the prevalent bi-conical structure 
of the AGN ionizati on in low-redshift active galaxies (e.g. 
iFischer et al.l[20T^ . The roughly vertical gap between 
the two blobs therefore marks the distribution of obscur¬ 
ing material with opacity higher than inside the [O ill] 
blobs. This picture is well consistent with the median 
velocity map on which the distribution of its highest 
and lowest velocities clearly depict a rotation pattern 
remarkably similar to a highly inclined {i > 70°) disk 
galaxy (to confirm this judgement, we use the publicly 
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Figure 1. IRAS F04250—5718. Top color maps. Maps of the [O iii]A 5007A surface brightness (left, in units of erg s ^ cm ^ arcsec 
logarithmic scale), median velocity (middle, in km s“^) and line width (VFso, right, in km s“^). The seeing at the observing site is depicted 
by the open circle on the surface brightness map. The white dashed line on the Rmed cind Wgo maps marks the orientation of the nearly 
edge-on galactic disk rotating at ~ 170 km s“^. Bottom spectra. We select 6 spatial positions to present the [O ill] velocity profile 
therein {Fx is in units of 10“^^ erg s“^ cm“^ within a 0.1" spatial pixel), including 2 outflow regions (“A”, “B”), 2 galactic disk 

regions (“C”, “D”), and two outer regions (“E”, “F”). The fitted line is in orange, the median velocity is marked by red dashed lines, and 
the velocity range used for calculating IFso is denoted by grey boxes. 


available software D i skFit (iSpekkens fc SellwQQdl2QQ7 ; 
iSellwQod fc SanchezI I2Q1QI: iKuzio de Narav et al.l 12012 ) 
to fit the velocity field of the central region to a rotating 
disk model, finding an inclination angle ofi = 70°it 10°). 
Half of the maximum velocity difference on the Vmed map 
is 170 km s“^, approximating the line-of-sight rotation 
speed of a typical edge-on disk galaxy, which extends 
along the direction ^ 20° from top to the right. 

The Wso maps reveal a pair of prominent high- 
linewidth lobes. This structure spatially coincides with 
the off-center brightness peaks seen in the [O ill] inten¬ 
sity map, and extends in a direction (nearly horizontal in 
the figure) ^ 60° from the direction of maximum Wmed 
difference which appears to fall in a low line-width valley 
(the contours marking approximately the maximum and 
minimum Umed values are overlaid on the Wgo map to 
guide the eyes). The line width of the pair of lobes reach 
as high as Wgo ^400-500 km s“^, strongly pointing to 
being highly tur bulent outflowing gas (cf. discussion in 
iLiu et ah l2Q13bD . 

Therefore, the three maps in Figure [1] are in consis¬ 
tency with a simple and familiar physical picture: the 
AGN resides in the center of a nearly edge-on disk galaxy 
with a rotation speed of ^ 200 km s“^, which both ob¬ 
scures its photoionization along the direction of the disk 


so that the surrounding gas forms a wide (> 120°) pho¬ 
toionization bi-cone, and forces the galactic wind to prop¬ 
agate into other directions with less obstruction. It is 
unclear why the outflow is tilted ^ 20° off the axis of the 
rotating disk. Possible explanations include a warped 
galaxy disk, or a detailed matter distribution with fur¬ 
ther complication that is not revealed by our observa¬ 
tions. 

High [O III] linewidth implies outflowing ionized gas, 
both because of the bulk motion of the gas, and more 
plausibly, because of the broad distribution of the veloc¬ 
ity of the nar row-line clouds em bedded in the winds (see 
discussion in iLiu et all l2Q13bD . In view of the uniform 
Hgo value inside the galactic disk (spanning a narrow 
range of 280-320 km s“^ over the narrow region where 
Vmed reaches its minimum and maximum. Fig. [T]), we 
conservatively choose the H80=320 km s“^ contour as 
representing the boundary of the galactic outflow. The 
furthest galactocentric distance that the outflow reaches 
is then measured to be 1.5", translating to 2.9 kpc, on 
the more extended side, and 1.1" or 2.1 kpc on the more 
compact side. 

These distances should be deemed as lower limits. The 
[O III] line emission is strongly detected in the entire 
10 X 10 kpc^ filed of view, thus the possibility cannot 
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be ruled out that Wgo falls below 320 km s“^ when 
the wind propagates to sufficiently large distances and 
even become undetectable. Further more, galactic winds 
are multi-phase and our observation only traces the 
warm ionized gas. A volume-filling, high-temperature, 
l ow density ionized g as phase has long been proposed 
(|Heckman et al.l llQQOh , which may be physically even 
more extended and cannot be traced by [O ill] emission. 
Projection effects also lead to underestimating the phys¬ 
ical scale if the outflow does not perfectly lie in the plane 
of the sky. In spite of these uncertainties, the basic con¬ 
clusion is validated that an outflow with a total physical 
extent of at least 5 kpc is exerting feedback effects on 
galaxy-wide scales. 

4.2. Mrk 509 

Mrk 509 is a more complicated system than F04250, 
and interpreting the [O ill] nonparametric measurements 
is nontrivial. In particular, t he HST FQ5Q8N narrow - 
band [O ill] imaging (Figure l^ lFischer et al.|[2QT^l2Q15l) 
detects tidal structures featured by a prominent linear 
tail extending out to a length of 7.8" (5.4 kpc) that 
appears to abruptly bend towar d the nucleus from the 
southwest. iFischer et al.l (|2Q13[) suggest that the posi¬ 
tive velocity of the southwest structure is indicative of 
inflowing, rendering the system a minor merger with a 
dwarf galaxy. These features are also spatially resolved 
by our IFU observations, which becomes more evident 
after careful PSF subtraction (detailed in Section [321) , 
and makes characterizing outflow features a challenging 
task (Figure |3l top row). 

The [O III] line shows unambiguous double-peak pro¬ 
files in the west ^40% fraction of the FoV. The two peaks 
are well separated wherever the tidal tail dominates the 
[O III] emission: the bluer peak always corresponds to a 
negative velocity (between —180 and —100 km s“^), and 
the redder peak persistently shows a positive velocity 
(70-220 km s“^). We therefore perform 2-Gaussian fits 
to the velocity profile of the [O ill] line in every spatial 
pixel of this region to decompose their respective spatial 
structures. 

This approach results in a clear and physical decom¬ 
position demonstrated in Figure |3l The intensity map 
of the Gaussian peaking at a bluer wavelength recon¬ 
structs the morphology of the linear tidal tail precisely, 
whose narrow linewidth (typically Wgo ^ 200 km s“^) 
throughout the region is also consistent with being an 
gaseous structure illumin ated by the centr al A GN (mid¬ 
dle ro w; see discussion in lLin et~ani2013bl and iLin et al.l 
I2009D . The intensity map of the re dder Gaussian r eveals 
a “southwestern jut” mentioned in IFischer et al.l (|2013D 
with little contamination from the tail (bottom row). 

Here we interpret the “southwestern jut” as part of the 
galactic outflow (see other possibilities and a discussion 
of the complex gas kinematics in Section 16.2p . This is 
not only because of the wholeness, smoothness and sym¬ 
metry of [O III] flux when it is mosaiced with the original 
data for the rest of the FoV (Figure [3l bottom row), but 
more importantly, because the ratio of [O ill] to the nar¬ 
row H/3 component is as high as 13 on average (same as 
that in the center) which indicates high ionization state 
caused by AGN emissi on routinely found in photoioniza¬ 
tion cones (see Section [6Tl iLin et al.ll2Q13aD . and because 
the mosaiced Wgo map shows behaviors similar to known 


outflowing objects (see below). 

Broad line widths are i ndicative of gaseo us outflows in 
IFU data, as reasoned bv lLin et al.l (|2Q13bD . in which the 
[O III] -emitting clouds possess a wide local velocity distri¬ 
bution in addition to the global outflow kinematics, lead¬ 
ing to broad wings of the [O ill] line profile. In absence of 
the tidal tail, we are left with a structure with elliptical 
or roughly round morphology shown by the [O III] flux, 
with relatively high Wgo values distributed in a round¬ 
ish region. The situation here is si milar to the lumin ous 
2 : ^ 0.5 SDSS quasars observed bv lLin et al.l (|2Q13bD . In 
those objects, high linewidths persist at a roughly con¬ 
stant level in the central several kpc (Wso ^ 500-2000 km 
s“^) which often decrease to much lower values (^100- 
200 km s“^) in the outer regions (several kpc from the 
center) beyond the direct impact of the outflow. 

In Mrk 509, Wgo declines rapidly at ^ 310 km s~^ in 
the eastern 60% of the FoV, and at ^ 220 km s“^ in the 
region where the tidal tail feature is removed through 
decomposition. Because 2-Gaussian fits cannot account 
for the broad base of the [O ill] line in this region, Wgo 
measured there may be underestimated. Whether or not 
taking into account this part of the FoV, we measure 
the northeastern extent of the outflow (where we are the 
most certain about its outflow characteristic due to the 
observed suppressed narrow H/3 line, see discussion in 
Section 0), and find a radius of 1.2 kpc (1.7") for this 
roughly spherical (or very wide angle) outflow. 

5. COMPARE TO ABSORPTION LINE ANALYSES 

5.1. IRAS FO 425 O- 57 I 8 

5.1.1. Spatial extent 

The galactocentric distance of the outflow revealed by 
our Gemini data, R ^ 2.9 kpc and R ^ 2.2 kpc for 
the two sides of the bi-conical outflow, is consistent with 
the previous UV absorption line analysis obtained by our 
group, R^3 kpc lEdmonds et ahl (|2QllD . In that paper, 
we measured the column density of C ll and an upper 
limit on the column density of C ll*, and determined the 
electron density in the outflow to be <30 cm“^. Gom- 
bined with photoionization modeling, we found the above 
constraint for the outflow radius. The direct IFU map¬ 
ping in this work provides the first ever direct test that 
validates our approach to indirectly derive the outflow 
size, which has bee n applied on nearl y 20 objects during 
the last decade (see lArav et al.|[2013l for a review). 

5.1.2. Gas Kinematies 

ILin et all (|2013bD proposed a spherical, extinction-free 
outflow model with a constant physical velocity as a rep¬ 
resentative picture for the quasi-spherical outflows seen 
in their luminous quasars at 2 : ^ 0.5. For such a sim¬ 
plistic model, the radial velocity profile remains exactly 
the same for the entire system (see Equation 10 therein), 
rendering the observed linewidth a constant across the 
Wso map. As discussed in Section m the outflow from 
IRAS F04250-5718 clearly shows bi-conical structures. 
The opening angle measured from the high lUgo regions 
(Figure H]) is ^ 70° for both cones. The spherical model 
is not directly applicable to this system, but the calcula¬ 
tion in that paper can be utilized here. 

As in Figure H] we do not observe an limb-brightened 
outflow (expected for a hollow-cone structure where 
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Figure 2. HST archival [O lll] narrow-band image of Mrk 509 (GO-12212; PI: D. Crenshaw). The yellow box depicts the 5"x5'' field 
of view of our Gemini IFU observation. The spectac ular 5.4-kpc-long linear “tail” structure and the “southwestern jut” that appears to 
connect it to the center are also marked. See Section IQ for discussion of these structures. 


the walls of the co nes dominate the emission; e.g. 
iCrenshaw et al.l[2QQQ[) . we assume a filled bicone model 
equivalent to part of a spherically symmetric outflow, in 
which the luminosity density only depends on the spher¬ 
ical radius r, and scales as a power law, j{r) oc r~^. The 
very high inclination of the galactic disk (see Section 0]) 
allows us to assume that the biconical outflow (roughly 
perpendicular to the galactic disk) to lie approximately 
in the plane of the sky. 

In this model, the maximum linewidth is achieved on 
the axis of the bicone. If we further assume a constant 
out flow velocity up, the calculation for a spherical model 
in (|Liu et al.ll2Q13bl Section 4.1, Equations 7, 8 and 9) 
remains valid. Although the conclusion therein that the 
radial velocity profile remains the same at every spatial 
position is no longer applicable, on the axis of the bi¬ 
cone, Equation (10) in that paper only requires a trivial 
modification. 


where / is the [O ill] flux at a given spatial position and 
a given wavelength corresponding to Vz , the line-of-sight 
component of the outflow velocity vq, R\\ is the projected 
2-dimensional radius from the center along the outflow 
axis (in the plane of the sky under our assumption), 0 is 
the opening angle of each cone (^ 70°), 

In the outer region (i^y ^1.3 kpc) of the quasar neb¬ 
ula along the outflow axis, the observed [O ill] surface 
brightness indeed closely follows a power-law relation 
/ oc therefore a = 4.0 as per our best fits. 

Inserting these parameters into Equation [H we And a 


Vo sin I 


(a-3)/2 


( 1 ) 


simple linear relationship between lEgo and uq, 

Wgo = 1.374 uo sin^. (2) 

The observed maximum Wgo value is found to be 498 
and 410 km s“^ on the axis of the larger and smaller 
cone, respectively. Due to the limitation of our spatial 
resolution, we estimate the possible range of the opening 
angle to be 60° ^ 0 < 90°, which introduces an ^20% 
uncertainty in the inferred uq. Eor 0 = 70°, we attain an 
outflow velocity of vq — 520-630 km s“^. 

If the observed Wgo values are predominantly due to 
the projected velocities of the bulk motion of the out¬ 
flow, the physical velocity of the outflowing gas must be 
high. The host galaxy, rotating at a speed of ^ 170 
km s“^ (measured at 1-2 kpc from the center, see Sec¬ 
tion 11]) close to that of the Milky Way (^220 km s~^ 
at th e corresponding galactic radius. iBhattachariee et ^ 
l2014l ). This rotation speed indicates a host galaxy 
with comparable or slightly smaller mass than the Milky 
Way (although the rotation we observe may not reflect 
the large-scale dynamical mass of the system, and this 
is only a heuristic comparison). Considering that the 
Milky Way has an escape velocity of 550.9122.1 
(jKafle et al.l l2nT4^ . the ionized gas traced by [O ill] is 
probably (barely) escaping from the potential well of its 
host galaxy. However, we emphasize that in addition 
to the bulk flow, the local velocity distribution of the 
clouds in the outflow likely also contribute to the ob¬ 
served linewidth, resulting in a smaller outflow velocity. 

Under the “Case B” assumption, the mass of the ob¬ 
served ionized gas is proportional to the product of H/S 
luminosity and electron density (jNesvadba et al.l I2Q11I: 
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Figure 3. Mrk 509. Top color maps: The [O ill] surface brightness (left column, in erg s ^ cm ^ A logarithmic scale except for the 
“tail”), median velocity (middle column, in km s“°^) and line width (Vkso, right column, in km s“°^) maps. Our 2-Gaussian decomposition 
of the [O III] line allows for isolating the contribution from the “tail” structure (Figure [2]) and a roughly spherical or elliptical structure 
likely due to an galactic outflow. Overlaid on the bottom right panel is the 220 km s“^ contour for the “tail” region and the 310 km 
s“^ contour for the rest of the FoV, which roughly depict the outflow structure. Red contour crudely marks where the narrow H/3 line is 
weakest (indicative of star formation suppression), which spatially coincides with the broadest [O ill] line (see discussion in Section [5] cf. 
Figure ID). Central AGN is marked with a black cross, and the seeing at the observing site is shown by open circles. Bottom spectra: [O 
III] velocity profiles at 6 representative spatial positions {Fx is in units of 10“^^ erg s“^ cm“^ A“^, within a 0.1" spatial pixel), including 
2 outflow regions (“A”, “B”), 2 outflow plus “tail” regions (“C”, “D”), and 2 outer regions (“E”, “F”). The fitted line is in orange, the 
median velocity is marked by a red dashed line, and the velocity range enclosing VFgo is a grey box. 
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lOsterbrock fc Ferlandl I2QQ6I: though this is a ctually a 
lower limit, see discussion in iLiu et alJl2Ql^ Sec. 6). 
In this object, H/3 follows the same spatial distribution 
as th at of [O ill] because [O lll]/H/3 is a constant (Sec. 
15.2.11 Fig. [5]). Assuming that the electron density is dis¬ 
tributed as Uq = 10^ (r/0.1kpc)“^ cm“^ (Sec. Ib.l.ip . we 
infer that ^60% of the warm ionized gas is outflowing, 
therefore likely escaping from the galaxy’s potential well. 

Absorption line analysis of the high quality HST COS 
spectrum reveal three kinematic components in this sys¬ 
tem, identified using the strong, unblended C IV and 
N V doublet lines. Their centroids are found at ra¬ 
dial ve locities of -38, -156 and —220 km s“^, respec¬ 
tively (jEdmonds et al.l[2Qll[ ). These absorption outflow 
components are all at velocities smaller than the max¬ 
imum line-of-sight velocity uq sin(^/2) of 300-360 km 
s“^ (Equation [2|) given by the bi-cone model built upon 
our lEU observations, implying crude consistency. The 
> 100 kms s“^ lower velocities in UV absorption are 
probably because these absorbers are far away from the 
primary biconical structure of the ionized gas outflow, 
indicating non-localized matter distribution and non- 
uniform velocity distributions beyond the scope of our 
simple model. These poorly understood uncertainties re¬ 
flect the limitations of current knowledge on the physics 
of outflows, more realistic modeling and further observa¬ 
tional scrutinization in the future are awaited to provide 
a handle. 

5.2. Mrk 509 
5.2.1. Spatial extent 

In order to compare the distribution of AGN activ¬ 
ity to that of newly formed stars, we decompose the 
H/3 line into broad (EWHM = 2700 km s“^) and nar¬ 
row (EWHM = 140-640 km s“^) components using 2- 
Gaussian fits, so that the emission from the broad line 
region and the narrow line region of the AGN is sep¬ 
arated (note that [O ill] is emitted by the narrow line 
region only). The maps of signal-to-noise ratios at the 
peak of [O ill] and narrow H/3 lines are shown in Eig- 
ure m where the tidal tail feature has been removed, 
as described in Section IH Eaint [O ill]-emitting gas il¬ 
luminated by massive stars are characterized by small 
linewidths (EWHM < 200 km s“^). Thus, the line in¬ 
tensity in those regions, scaling as the product of the 
linewidth and the peak flux, are thus more suppressed 
in [O III] flux maps than in these signal-to-noise ratio 
maps. As a result, these maps are better for scrutinizing 
low surface brightness features than line intensity maps. 

In figure m at least 5 [O ill]-emitting blobs surrounding 
the center are seen, which are distributed at a galacto- 
centric distance of about 1.6 kpc. All these features are 
clearly seen in the narrow H/3 signal-to-noise map at the 
same positions. In galactic environments, hydrogen re¬ 
combination lines are emitted by interstellar g as ionized 
by you ng, massive stars with ages <10 Myr (|Liu et al.l 
l2013c[ ). rendering th em reliable trace r s of the current 
star formation rate (|Kennicu t3 119981: iLiu et all l2011f) 
if dust extinction is carefully corr ected for spatially- 
reso lved studies (see discussion in iGalzetti et al.l 120071 
and ILiu et al.ll2013dh . Therefore, this spatial correspon¬ 
dence of [O III] and narrow H/3 implies significant contri¬ 
bution or even dominance of massive stars in producing 


the warm ionized gas at this distance. 

Eurthermore, the ratio of these two maps reveal that [O 
lll]/H/3narrow IS abovc 10 in the central ^ 1 kpc and de¬ 
creases to 1-3 at a radius of 1.6 kpc and b eyond (see Sec- 
tionld.ll and Eigure[5]). In a BPT diagram ([Baldwin et al.l 
119811) . [O IIl]/H/3narrow > 10 implies AGN dominance, 
while [O lll]/II/3narrow ^ 1 points to a low-ionization nu¬ 
clear emission-line region (LINER, not applicable to this 
case) or star formation possibly mixed with AGN ioniza¬ 
tion. The combination of this fact and the spatial cor¬ 
respondence of [O III] and narrow H/3 features strongly 
indicate that the AGN’s dominating role in the interstel¬ 
lar radiation is taken over by star formation at a distance 
of ^ 1-1.5 kpc. We therefore conclude that the full ex¬ 
tent of predominating AGN ionization likely have been 
probed, in spite of our limited field of view (3.3 x 3.3 
kpc^). 

An intriguing phenomenon is the indicative evidence 
for the outflow suppressing star formation: the weakest 
H/3 (narrow) emission is found in an elongated region in 
the northeast to north direction from the center, at a ra¬ 
dius a ^ 1 kpc (delineated with a cyan dashed curve in 
Eigure[5]); this region spatially coincides with where the 
line width Wgo reaches its maximum (Eigure[3]). Similar 
evidence for outflows suppressing star formation activi¬ 
ties has be en reported in spatially resolved spectroscopy 
at ^ = 2.4 (|GanQ-Dfaz et al.ll2Q12| ) . This further consoli¬ 
dates the conclusion in ILiu et al.l (|2Q13bD that high line 
widths of [O III] strongly implies outflowing gas. 

The multi-wavelength campaign on Mrk 509 has de¬ 
rived a lower limi t of 100-200 pc for the outflow ra¬ 
dius in this system lArav et aT] (|2012D . This was inferred 
from the fact that the column densities of C IV and N v 
showed negligible variation between 2001 and 2009 de¬ 
spite a large change in ionizing flux, aided by Monte 
Garlo simulated light curves that statistically determined 
the distance limits. The galactocentric distance of 1.2 
kpc that we determine from the features of a spherical 
outflow (Section [4]) lies above these indirectly inferred 
lower limits comfortably. 

5.2.2. Gas Kinematies 

The ex tinction-free quasi- spherical outflow model pro¬ 
posed in ILiu et all (|2Q13b[ ) is applicable to Mrk 509, 
because of the round-ish morphology, the small veloc¬ 
ity variation across the median velocity map, and the 
roughly flat radial profile of Wgo within a radius of ^ 1 
kpc. Similar to the case of IRAS E04250-5718, we mea¬ 
sure the outskirt of the ionized gas nebula (before reach¬ 
ing the star formation dominated region at R > 1.3 kpc), 
finding the best-fit power law of [O ill] surface brightness 
to be / oc Ghanging Equation [1] to its orig¬ 

inal form in ILiu et ahl (|2Q13b[ ) (i.e. replacing Ry and 
vq sin(^/2) with R and uq, respectively), we find for 
a = 3.9, 

IFgo = 1.393 uq. 

To avoid additional uncertainties introduced by removing 
the linear tail structure, we only use the west (left) 60% 
part of the Wgo map (Eigure [3|) where nonparametric 
measurements were performed on the original [O ill] line 
profile. The average Wgo of the outflow in this region 
(defined by Wgo > 310 km s“^) is 410 ± 72 km s“^. 
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Figure 4. Maps of signal-to-noise ratios measured at the peak of [O lll] 5007A and the narrow H/3 emission line profiles. These maps 
allow for seeing low surface brightness features in the outer region more easily than previous flux maps (Figure O left panels). At least 
5 [O lll]-emitting blobs surrounding the center are seen in both panels (marked “E” to “I”). Note that the spectra at “E” and “F” have 
been shown in Figure [3] The cyan contour depicts the region of the weakest narrow H/3 emission, which spatially coincides with where VFso 
reaches its maximum values, indicative of star formation being suppressed by the outflow (cf. Figure [3] Wgo maps). 


translating to a physical outflow velocity of vq = 293 ±51 
km s“^, where the standard deviations are reported as 
uncertainties. Wgo reaches it maximum ^ 510 km s“^ 
in the northeast of the AGN, thus the largest possible 
outflow velocity is vq ^ 370 km s“^. 

Since this outflow is quasi-spherical, velocities from ab¬ 
sorption line analyses can be directly compared to our 
model (projection effects are minimized). Nine indepen¬ 
dent outflow components have been identified fro m the 
N V, Si IV and C IV troughs bv lArav et al.l (|2012[) . The 
strongest absorption feature emerges at r\j — 320 km s 
in agreement with our best estimate vq = 293 ±51 km 
s“^. The highest velocity measured in that work (425 km 
s“^) slightly exceeds the upper bound we derive here (370 
km s“^), but the corresponding component is relatively 
weak. Moreover, as discussed in Section 15.1.21 absorp¬ 
tion line analyses intrinsically have significantly higher 
spatial resolution and are more sensitive to the local ve¬ 
locity distribution of the clouds residing on our line of 
sight, which is likely smeared out in the IFU data which 
may only detect the bulk flow. Besides the uncertainty of 
our simplistic model, this difference could also be respon¬ 
sible for the small discrepancy (several tens of km s“^). 
As a result, these two independent and complementary 
investigations provide consistent results. 

5.3. Electron density 

Our IFU spectra include the conventional electron den¬ 
sity (ue) diagnositic, [S ll] AA6717, 673lA, facilitating a 
comparison of Uq to our previous absorption line analy¬ 
ses. To do this, we fit the [S ll] doublet simultaneously 
using the same line width and velocity structure to derive 
a [S ll] 6731A/6717A ratio map. 

The theoretical dependence of this line ratio on is 
then calculated using Version 7.1 of CHIANTI, an atomic 
database for spectroscop i c diagnostics of a strophysical 
plasmas (|Dere et al.lll997l:rLandi et al.ll2Q13[) . In this cal¬ 
culation, we assume a typical temperature of warm ion¬ 
ized gaseous nebulae that produce narrow emission lines, 
T = 10^ K (but note is only weakly dependent on the 


assum ed temperature as oc T ^/^. lOsterbrock fc Ferlandl 

[200l) . 

As a result, we find = 1300 ± 140 cm“^ and 
1200 ±110 cm“^ in the center of IRAS F04250—5718 and 
Mrk 509, respectively. We limit our measurement to 
locations where the [S ll] doublet has a signal-to-noise ra¬ 
tio of at least 3. At the maximum radii set by the above 
requirement, we find < 100 cm“^ at 1 kpc away from 
the center of Mrk 509, and < 200 cm“^ at 2 kpc from 
the center of IRAS F04250—5718 (these values are up¬ 
per limits because already reaches the lower critical 
density for deexcitation, where the line ratio is no longer 
sensitive to Ug). In Section EH Ue at the break radius 
is necessary for deriving the energetics of these outflows, 
which is, unfortunately, not measurable from our data, 
but the adopted value of 10 cm“^ is consistent with these 
limits. 

In the absorption outflow of IRAS F04250—5718, a 3cr 
upper limit of < 30 cm“^ is derived from the con¬ 
straint on the C II* /C ll ra t io pla ced by the non-detection 
of C II* (|Edmonds et al.l l2011[) . These absorbers are 
found to be located at ^ 3 kpc from the center; at this 
distance, we derive a consistent upper limit from [S ll] 
emission lines, < 200 cm“^^ _ 

UV absorption lines allowed lArav et ^ (|2012[) to mea¬ 
sure electron density for outflow components T1 and T2 
in Mrk 509, and cm“^, respec¬ 

tively, but the distance for the absorbers is only con¬ 
strained to be > 200 pc. At these distances, the [S ll] 
emission lines in our data giv e rig < 900 cm s~ ^, a re¬ 
sult in broad agreement with lArav et ^ (|2012[) , given 
the ^ 10% uncertainty. 

6. DISCUSSION 
6.1. Outflow energetics 

In a possibly realistic physical picture, numerous rela¬ 
tively dense clouds filling only a fraction of the volume 
as low as < 1%, which are embedded in a rarefied hot 
wind which might even dominate the mass and energy 
of galactic outflows. At small galactic radii, the sur- 
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face of these clouds are ionized by the central AGN and 
emits the observed narrow lines (“ionization-bounded”). 
As the wind carries these clouds to larger radii, the de¬ 
creasing pressure allow them to expand and gradually 
become optically-thin, and finally fully ionized at a cer¬ 
tain distance , entering the “ni atter-bounded regime” (cf. 
discussion in lLiu et al.ll2Q13'bh . 

This hypothesis of transiti on is supported by a 
series of recent ob s ervat ions (jLiu et al.l l2Q13ai l2Ql4 
iHainline et al.l I2Q13LI2Q14D , of which the most direct is 
the radial profile of [O lll]-to-H/3narrow persisting at > 10 
at small radii that starts to decline at a break radius 
of Rhr = 4-11 k pc, beyond which He n/R/3 increases 
(iLiu et al.ll2Q13aD . This newly found phenomenon has 
enabled iLiu et all (|2Q13b[) to constrain the estimation of 
the mass flow rate (M) and kinetic luminosity (^kin) of 
the outflows significantly better than previous studies. 
In this section, we follow that work to constrain the en¬ 
ergetics of the two objects. 


6.1.1. Energetics: Mrk 509 

The above [O lll]/H/3 behavior is also seen in Mrk 509, 
where the break radius is i?br — 1.1 kpc (Figure [5j). At 
this radius, we measure the H/3 surface brightness to be 
^H/3 = 2.6 X 10“^^ erg s“^ cm“^ arcsec“^ after correcting 
the cosmological dimming effect by a factor of (1 + z)^. 
This obse rvation fa c ilitates direct application of the rea¬ 
soning in ILiu et all (|2Q13b[ ). 

The electron density in outflows is poorly known. 
In mo re luminous and en ergetic objects at higher red- 
shifts, |Greene^F^ (|2^3)^ndvM00_cm“^_m radio- quiet 
quasar nebulae, while iNesvadba et al.l (j2QQ6l . l2QQ8l ) find 
a few 100 cm“^ in powerful radio galaxies. These results 
are not dissimilar to the densit y in the sub-kpc neigh¬ 
borhood of our AGNs fSection IO|) . The observed [O 
III] surface brightness profile, scaling as R~^ in both our 
objects (Section 15.1.21 an d Section 15.2.2]) , is a natural 
consequence of oc r~ ^ (|Liu et al.ll2013a[ ). In both ob¬ 
jects, the [S ll] emission lines give ^ 1000 cm“^ at a 
galactocentric radius of 0.1 kpc; it will fall to ^ 10 cm“^ 
at the break radius (1.1 kpc for Mrk 509). Although we 
can only measure an upper limit at ^ 1 kpc {ue <100 
cm“^; Section 15.31). this e xtrap olation leads to sim¬ 
ilar to iRupke fc Veilleu:?d (|2013D : those authors assume 
rie = 10 cm“^ based on measuremen ts of the supperbub- 
ble in NGC 3079 (rze = 5-100 cm~^: r\Ailleux et fflll994 
iGecil et ahl 120011) . The above considerations lead us to 
adopt Ue = 10 cm“^ for the following calculations. 

We find the line emission to fall off as I;h /3 oc 
in this object, changing the n umerical factor 1.44 in 
Equation 18 in ILiu et al.ll2013bl to 1.94. Using Equation 
22 therein, we estimate the mass flow rate to be 


M 


^YL^iRhT 


5.2MQyr ^ 2.6x10 ^^ergs ^ cm ^ arcsec ^ 


Rhr 


Vo 


1.1 kpc 293kms“^ VlOcm 




which is normalized by our derived outflow velocity, 293 
km s“^ (Section [5]). As a consequence, the kinetic lumi¬ 


nosity of the outflow, Ekin = Mvq/2^ is given by 

-^kin ^ M / Vq 

1.4 X 10^^ ergs“^ 5.2MQyr“^ \293kms“^/ 

The momentum flow rate isP = Mu = 9.6x 10^^ dyne, 
or log (cP/Lq) = 10.9, similar to th at of the outflows 
driven by low-redshift major mergers (|Rupke fc Veilleul3 
(|2013[ ) find a range of 10.0-12.8). 

The uncertainty of these results is largely caused by 
the poorly constrained ng. M and Pkin for other electron 
densities can be easily calculated using the n~^ depen¬ 
dence (Equations [3] and HI). The conservative constraint 
that [O III] -emitting clouds only fill a fraction of the vol¬ 
ume results in a lower limit of Uq ^ 0.1 cm“^ (Equation 
21 in ILiu et al.ll2013'bl) . Along with the measured upper 
limit, Uq < 100 cm“^ fSection 15.3p . we find firm lower 
limits using Equations [31 and IH 

0.5 Mq yr“^ < M < 520 Mq yr“^, 

1.4 X 10"^^ erg s”^ < Pkin < 1-4 x 10^^ erg s-^ 

Hence, we conclude that Mrk 509 is driving an outflow 
with M ^ 5 Mq yr“^ and E ^ 1 x 10^^ erg s“^, with 
an upper (lower) limit of 100 (10) times larger (smaller) 
for both quantities. The ratio of the momentum flux of 
the outflow to the AGN radiation is P/(I/boi/c) ^ 0.3. 
The outflow kinetic luminosity is only ^ 0.01% of the 
bolometric luminosity of Mrk 509, and ^ 0.002% of 
its Eddington lumi nosity (Phni = 1 x 10^^ erg s“^, 
Pboi/AEdd — 0.16; iKaspi et al.l l2000h . Since theoreti¬ 
cal modeling predicts that significant feedback requires 
> 0.5-5% of the AGNs luminosity to be converted to the 
mechanical energy of the ou tflow (e.g. iDi Matteo et al.l 
l2005l : iHopkins fc ElvisI 120101 ). the outflow in Mrk 509 
does not produce feedback strongly affecting the evolu¬ 
tion of the host galaxy. 

On the absorption line front, lArav et al ] (HoHI) find 
an average hydrogen column density of A^h = 10^^’^ 
cm“^ and a velocity of ^ 300 km s“^ in the outflow 
of Mrk 509. Adopting a solid angle of Air for the quasi¬ 
sphere geometry and an outflow radius of 1.2 kpc ob¬ 
tained in this work, we find M = 3.3 Mq yr“^ (Equation 
T. lEdmonds et al.ll2011l ). and Pkin = Mu^/2 = 9.3 x 10^^ 
erg s“^, in agreement with the above determination 
within a factor of 1.6. 

6.1.2. Energetics: IRAS E04250—5718 

The uniform [O lll]-to-H/3narrow ratio > 10 across 
the whole map implies a break radius beyond the 
largest radius reached by our field of view (Pbr ^ 5.5 
kpc, see Eigure However, the largest break radius 
that has ever been observed is 11 k pc, seen in SPS S 
J032144.11+001638.2 at z = 0.643 (It4ii et, a, 1.1 [201^ . 
one of the most luminous quasars at that redshift, with 
a bolometric luminosity of > 10^^'^ erg s“^. IRAS 
E04250—5718 is expected to have Pbr ^ 10 kpc due to its 
lower luminosity (Pboi ^ 9x 10^^ erg s~^. lEdmonds et~al] 

Imll). 

At the 5.5 kpc radius, the [O ill] surface brightness 
is measured to be 1.3 x 10“^^ erg s“^ cm“^ arcsec“^ 
after correcting for the cosmological dimming effect, thus 
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Figure 5. The radial profile of the intensity ratio of [O ill] 5007 to the narrow component of H/3 (obtained from 2-Gaussian decomposition) 
for IRAS F04250—5718 and Mrk 509. The ratio persists at a constant level (~10) in the whole field of view for the former, bu t starts to 
declin e rapidly at a “break” radius of 1.1 kpc for the latter (a prevalent phenomenon discovered in 2 ; ~ 0.5 luminous quasars bv ILiu et al.l 
I2013alb This radius (marked with a vertical line) implies a transition of the gaseous nebulae from an ionization-bounded regime to a 
matter-bounded regime, and facilitates our calculation of the kinetic luminosity and the mass flow rate of the outflows. 


the surface brightness profile Il[oiii] oc i^-3-0±o.i (^Section 
I5.1.2P predicts S[oiii](ll kpc) = 1.7 x 10“^^ erg s“^ cm“^ 
arcsec“^, and thus i;H/ 3 (llkpc) > 0.1 H[oiii] = 1-7 x 
10“^^ erg s“^ cm“^ arcsec“^. 

Compared to a spherically symmetric geometry, the bi- 
conical configuration of the outflow introduces a correc¬ 
tion factor of 1 — cos (^/2), where ^ 70° is the opening 

angle of each cone (Section |4]). Now that both S[oiii] and 
I;h /3 fall off as the num erical factor 1.44 in Equa¬ 
tion (18) in iLiu et"^] (j2013b[ ) increases to 1.57. For an 
outflow velocity of 520 km s“^ (the lower bound, Sec¬ 
tion [51121), assuming i?br = 11 kpc and Ue = 10 cm“^. 
Equation [3] gives 

M > 1.3 Mq yr“^, 

P > 4.3 X 10^^ dyne, Pkin > 1-2 x lO"^^ erg s“\ 

As discussed in Section (6.1. 11 compared to Mrk 509, this 
object has similar Uq in the center and similar Uq radial 
profile, but its break radius is larger by a factor of a few 
or more, and at the break radius is likely lower than 
10 cm“^. However, only an upper limit, Uq < 200 cm 
s“^, can be derived from our data, after all. Here we 
adopt He = 10 cm“^ for consistency (which pushes the 
above results to their lower limits, see below). 

Since the actual electron density is likely lower, and the 
smallest possible Eh /3 is used here; although Pbr may 
have been overestimate by a factor of ^ 2, the overly 
large break radius leads to an underestimate of the steep 
declining Eh /3 (by a factor of ^ 8) that itself cannot 
balance out. For the above reasons. Eg nation 16.1.21 gives 
lower limits of the outflow energetics. 

On the other hand, adopting Uq = 0.1 cm“^, Pbr = 5.5 
kpc and vq = 630 km s“^ (upper bound) sets conserva¬ 
tive upper limits 

M < 690 Mq yr“^, 

P <2.7 X 10^^ dyne, Pkin < 8.6 x lO"^^ erg s^^ 

The ratio of the momentum flux of the outflow 
to the AGN radiation is only loosely contrained to 


P/(I/boi/c) = 0.01-9. The kinetic luminosity of the 
outflow in this system is insignificant, being only > 
0.002% of its bolomet ric luminosity (^ 9 x 10^^ erg s“^, 
[Edmonds et al.ll2QTTI) . Therefore, the outflow in IRAS 
F04250—5718 is not a significant feedback agent, similar 
tc ) the case of Mrk 509. 

[Edmonds et al.[ (|2011[) And a velocity of 220 km s“^ 
and a column density lower limit of A^h > 10^^'^^ cm“^ 
in the absorption outflow of this object. Adopting an 
outflow radius of 3 kpc, a solid angle of O/Att = 0.2 de¬ 
rived from our IFU data and deprojecting their velocities 
using 0 = 70° (the outflow lies approximately in the plane 
of the sky. Section [6.1.2p . their results lead to M ^ 1.2 
Mq yr“^, and Pkin ^ 5.5 x 10^^ erg s“^, in agreement 
with the above determination within a factor of ^ 2. 

6.2. More on Mrk 509’s structure 

Recent long-slit ([Fischer et al.[ [2013[) and IFU spec¬ 
troscopy (|Fischer et al.[[2015[) demonstrate that Mrk 509 
is likely a minor merger system, where gas is inflowing 
and fueling the central AGN through the linear tidal tail 
and the “southwestern jut” (cf. Figure [2]). In this sec¬ 
tion, we discuss in detail the dynamical structure of this 
system by combining and comparing our data to previous 
observations. 

The host galaxy is nearly face-on {i = 36?4 as per the 
Hyper Leda databas^l, [Paturel et al.[[200^ [Yaaoob et al.[ 
[2QQ3[ And i = 41°), and the linear tidal tail is likely lo¬ 
cated in front of the disk, as otherwise it would have been 
severely extincted by dust and become invisible. (To test 
this judgement, we isolate the tidal tail’s emission in the 
Ra and R/3 lines using 3-Gaussian fits accounting for the 
broad base and the double-peak structure, and And the 
average Ha/H/? ratio on the tail to be 2.8 ± 0.5. Un¬ 
der “Gase B” assumption, the intrinsic line ratio is 2.7- 
3.0 for a wide range of temperature and electron density 
(T = 5000-30000 K, = 10^-10^ cm-^). Hence, the 
observed line ratio is consistent with negligible dust ex¬ 
tinction, confirming our conclusion that the galactic disk 

^ http://leda.univ-lyonl.fr 
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is behind the tidal tail.) Therefore, in case the “south¬ 
western jut” is physically connected to the tail and serves 
as a tunnel for inflowing gas, we expect to see velocity 
smoothly transition from ^ —100 km s“^ to zero then to 
positive velocities, along the route of the tail, then the 
“jut”, finally the AGN neighborhood. 

Nevertheless, in our data we find that the “jut” and the 
long linear tail, although appear to be connected, have 
distinct velocities everywhere (different by 200-350 km 
s“^). Furthermore, the [O ill] intensity map of the tidal 
tail (Figure [3]) , middle row) shows that it likely switches 
direction sharply (southward to eastward) to approach 
the nucleus before it reaches the main body of the “jut”, 
at only ^ 0.4" (R.A. distance) south of the AGN. This 
phenomenon i s confirmed by IFU maps with higher spa¬ 
tial resolution iFischer et al.l (j2015f ). Hence, we conclude 
that the linear tail and the “jut” are physically unrelated, 
but are superposed by projection e ffects; the gas in th e 
linear tail is likely indeed inflowing (jFischer et al.ll20l'^ . 
but it may use its own route instead of inflowing through 
the “southwestern jut”. 

However, with the linear tidal tail removed using 2- 
Gaussian fits, interpreting the rem aining structures (in- 
cluding the “jut”) is still nontrivial. IFischer et al.l (j2015l ) 
performed IFU spectroscopy assisted by adaptive optics 
in the [S ill] 0.95 /im emission line. In their data with 
higher (^0.1") angular resolution, they observe a ^2" 
elliptical region (schematically shown by a ellipse in Fig¬ 
ure El in their smaller FoV (3"x3"). This region shows a 
well defined velocity gradient (upper left panel of Figure 
4 in their paper), allowing them to fit the velocity field 
to a rotating disk model. 

However, this rotation pattern is not seen in our me¬ 
dian velocity maps (Figure [3]): we do observe negative 
velocities close to the east edge of the FoV, but we see a 
roughly constant velocity across a large round-ish region 
when the tail is removed, which is deemed as indicative 
evidence for a quasi-spherical outflow in Section 14.21 

This seeming discrepancy is, in fact, due to the higher 
signal-to-noise ratios of our data. In fact, if we determine 
the velocity field from the peak wavelength of the [O 
III] line instead of calculating the m edian velocity, the 
pattern seen by IFischer et al.l (j2Q15[ ) is reproduced with 
accuracy. As shown in Figure [6l the rotation axis of their 
best-fit model (dashed line) is consistent with our zero- 
velocity spatial pixels (here we force the velocity at the 
center to be zero for comparison purposes). Obviously 
the rotation pattern is actually not restri cted within a 
circumnuclear disk (^0.7 kpc in radius) as IFischer et al.l 
(l2nTa find, but exists on a scale as wide as 3.3 x 3.3 kpc^ 
(our whole FoV) or larger. This finding may alleviate 
the issue of the rather high rotation speed those authors 
derived (^ 500 km s“^), if their model fits are performed 
on a n area larger by a fa ctor of ^ 5. 

In iFischer et ahl (f2015f ). [S III] is detected with S/N < 3 
in ^ 40% area of the FoV (roughly between the ellipse 
and the box in Figure [6]), where we detect [O lll]5007A 
line with S/N > 15 in every spaxel (Figure |4]). In low 
S/N data, the weak broad base of an emission line, which 
may contain outflow signatures, can be easily submerged 
into noise, so that a single Gaussian is sufficient for line 
fitting (this effec t has been simulated and discussed in 
I Liu et ahl (|2013bl . Section 3.4). This difference explains 



Figure 6. [O lll]5007A emission line’s peak velocity map of Mrk 
509 (in km s“^; the velocity at the AGN position is set to be zero). 
In contrast to the median velocity map (more sensitive to the broad 
base of the line profile, Figure [3|) that shows a roundish structure 
in the center likely caused by an outfl ow, this map exhibit s a well 
defined large-scale rotation pattern. IFischer et ahl (|2015l i And a 
similar rotating pattern in the [S ill] 0.95 p,m emission line, with 
an axis (dashed line) in consistency with this map; those authors 
attribute this motion to a rotating disk ~ 2" in diameter (schemat¬ 
ically shown as an ellipse) in the center of their 3" x 3'' FoV (box), 
but we see the rotation taking place on a larger scale, co-existing 
with a quasi-spherical or wide-angle outflow. However, the uncer¬ 
tainties that the rotating gas intr oduces into our outflow modeling 
are insignificant (see Section [6. 21) . 

why in majority of their spaxels [S III] can be fitted with 
a single Gaussian profile, which only happens in a few 
percent of the spaxels in our [O ill] analysis. It also re- 
sul ts in our veloc i ty ma p different from the one obtained 
bv IFischer et al.l (j2Q15f ). as median velocity is more sen¬ 
sitive to the broad base than fits to lower S/N data, and 
thus likely tra ces outflow ki n emati cs better. 

Meanwhile, IFischer et al.l (j2Q15[ ) also detect a second 
Gaussian component northeast to the AGN, which they 
suspect is an inflow or outflow. We now obtain sup¬ 
portive evidence for an outflow — this region spatially 
coincides with where [O ill] is broadest and where the 
narrow H/3 is weakest (Figures [3] and [5]), indicating star 
formation being suppressed by an outflow. This north¬ 
east outflo w region show s an o pening angle of ^ 90° in 
the data of IFischer et al.l (j2Q15l ) , whereas our data points 
to a significantly wider angle or even a spherical morphol¬ 
ogy (Section l4.2[) . 

The above considerations lead us to conclude that a 
wide-angle outflow and a large-scale rotating gas likely 
both contribute to the kinematics (beside the inflowing 
tidal tail in the west). Decomposing the two components 
is difficult, either spatially or spectroscopically, but we 
reckon that the contribution from the rotating gas is not 
significant. Empirically, rotating gas in the host galaxy 
produces emission lines narrower than those from out¬ 
flows (HVo ^ 200 km s“^ vs. HVo ^ 300-400 km s“^). 
Therefore, in general, Hgo directly measured from the 
observed line profile is not significantly different from 
the line width of the outflow only (unless the outflow is 
remarkably fainter than the rotating gas, in which case 
the measured Wgo provides a lower limit). This as s ertion 
is supported by the existing data: IFischer et al.l (j2015l ) 
separate out the outflow contribution using two-Gaussian 
fits and finding the outflow component in the northeast 
region to have HVo = 500-600 km s“^, in agreement 
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with our direct measurement from the observed [O III] 
line (Figure [3]). Hence, the line width used for calcu¬ 
lating the outflow velocity in Section 14.21 is confirmed by 
an independent observation, and the contamination from 
the rotating galactic gas is trivial. 

Thus far, we have assumed a spherical geometry for 
the outflow in Mrk 509. As discussed above, the north¬ 
east quadrant is where the existence of an outflow is 
supported by the strongest evidence. Since the firm 
l ower limit of its o pening angle is ^ 90° as given by 
([Fischer et al.|[2QT5l) . assuming such a biconical geome¬ 
try with its axis close to the plane of the sky leads to a 
velocity larger by ^ \f2 and an additional filling factor of 
~ 0.3, and thus a mass flow rate and a kinetic luminosity 
not dissimilar (within a factor of ^ 2) to the results of 
assuming a spherical symmetry. 

7. SUMMARY 

This work is the first effort to bridge integral field spec¬ 
troscopy and absorption line analysis investigations on 
the ionized phase of AGN outflows, the two indepen¬ 
dent and complementary approaches to determine out¬ 
flow properties. We present Gemini IFU observations 
of the ionized gas nebulae around two low-luminosity 
quasars, Mrk 509 and IRAS F04250—5718. The outflow 
distances from the central AGNs have been indirectly 
determined from UV absorption line analyses, which are 
directly tested by our IFU observations in this work. 

We demonstrate that IRAS F04250—5718 shows a clear 
case where a biconical outflow extending out to at least 
2.9 kpc from the center, in agreement with the result of 
absorption line analysis (R ^ 3 kpc). The opening angle 
of each cone is ^ 70°, and the outflow axis is ^ 20° away 
from the axis of the galactic disk rotating at a speed 
of ^ 170 km s“^. The outflowing gas, with a derived 
physical velocity of 580 ±80 km s“^, is probably (barely) 
escaping from the potential well of its host galaxy. We 
obtain lower limits of the mass flow rate, M > 1.3 Mq 
yr“^, the momentum flux rate P > 4.3 x 10^^ dyne, and 
the kinetic luminosity ^kin > 1.2 X 10^^ erg s 

Mrk 509 is a more complicated system especially be¬ 
cause of a linear tidal tail. After removing the tidal tail 
spectroscopically, we find a quasi-spherical outflow with 
a radius of ^ 1.2 kpc, also consistent with pre-existing 
absorption line analysis (a lower limit of 100-200 pc). 
The spatial coincidence of the highest [O ill] line widths 
and the weakest narrow H/3 emission is indicative evi¬ 
dence for the outflow suppressing the formation of stars. 
Using the derived outflow velocity of 290 ± 50 km s“^, 
we estimate a mass flow rate of M ^ 5.2 Mq yr“^, a mo¬ 
mentum flux rate of P ^ 9.6 x 10^^ dyne, and a kinetic 
luminosity of Pkin ^ 1.4 x 10^^ erg s“^. 

Adopting the outflow radii and geometric parameters 
measured from IFU, absorption line analyses would yield 
mass flow rates and kinetic luminosities that are in agree¬ 
ment with our IFU determination within a factor of ^ 2. 
Measuring the outflow properties is challenging, and the 
multiple difficulties lead to relatively large uncertainties. 
However, we emphasize that the lower and upper lim¬ 
its we derived in this work spanning wide ranges are 
rather (maybe even overly) conservative, and should be 
understood as firm limits instead of limits with statisti¬ 
cal meanings (e.g. 3cr limits). The measured values of 


the outflow rates and energetics are reasonably reliable, 
and the above agreement between IFU and absorption 
line analyses is a justified validation. 

The [S ll] AA6717,6730 A doublet allows us to derive 
electron densities in the central 1-2 kpc of both objects. 
The derived electron densities are in agreement with the 
previously inferred values from their absorption outflows 
(Mrk 509: Uq ^ 10^ cm“^ at ^100 pc from the center; 
IRAS F04250-5718: < 200 cm-^ at > 2 kpc). 

We conclude that both the spatial locations and the 
kinematics of these outflows determined from IFU spec¬ 
troscopy of emission lines are consistent with those de¬ 
rived from previous UV absorption line analyses. The¬ 
oretical modeling predicts that significant feedback re¬ 
quires ^ 0.5-5% of the AGN’s Eddington luminosity 
to be converted to the mechanical energy of the out¬ 
flow. Feedback in these systems is taking place on kpc 
or galaxy-wide scales, but is inadequate for regulating 
the evolution of their galaxy hosts effectively (kinetic lu¬ 
minosity is only ± 0.002-0.01% of the bolometric lumi¬ 
nosity), unless it used to be significantly more energetic 
in the past. 
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